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Non-linear CMB temperature anisotropies up to the third-order on large scales are calculated. On large 
scales and in the Sachs- Wolfe limit, we give the explicit expression for the observed temperature anisotropy in 
terms of the primordial curvature perturbation up to third-order. We derived the final bispectrum and trispec- 
trum of anisotropies and the corresponding non-linear parameters, in which the contributions to the observed 
non-Gaussianity from primordial perturbations and from the non-linear mapping from primordial curvature per- 
turbation to the temperature anisotropy are transparently separated. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 



In the past few years, extensive attention has been attracted to the investigation of cosmological perturbations yD beyond 

the Unear order. The importance of studying non-linear perturbations comes two aspects. Firstly, forthcoming experiments 

^^ . of Cosmic Microwave Background (CMB) and Large-scale Structure (LSS) will be able to detect the non-linear structures in 

f) these perturbations. The observational detection of the non-linearities through the statistical non-Gaussianity (see 14J] for a 

comprehensive review and 121 for a review of recent observational progress ) of perturbations has become one the the primary 

targets of the cosmology. On the other hand, non-linearities, which encode the interactions in the early universe, would definitely 

bring us new understandings of both the early universe and the fundamental physics. 

O ■ A large amount of efforts have been devoted to the calculation of the statistics of curvature perturbation (, like primordial bis- 

|H ^ pectrum and trispectrum on large scales, pioneered by Maldacena 115711 ' . However, these are not the observed non-Gaussianities 

^ i of (e.g.) CMB temperature anisotropies ^ . Conventionally, one may use the linear order relation -^ = — ^$ = — |Cto eval- 

i ' uate the angular bispectrum or trispectrum of CMB, assuming the contributions from the second-order or secondary effects are 

negligible comparing to the primordial ones llll2,lal- However, in light of increasingly precise observations, a full treatment of 

the higher-order radiation transfer functions of the CMB anisotropies is needed, which will allow us to make definite prediction 

^^ ' of CMB non-Gaussianities (see 15|l for a recent review for non-Gaussianity on the CMB). 

^ ' The research on non-linear temperature anisotropy due to gravitational perturbations was pioneered by lIs UlOll . in which the 
^vj , second-order generalization of Sachs-Wolfe (SW) effect and Integrated Sachs-Wolfe (ISW) effect were derived. These results 
,__! . were extended in lll3l [l4l [l8l - l20ll where the second-order radiation transfer function on large scales was calculated, and in lllSll 
• ■ where the general expression for anisotropy due to gravitational perturbations up to the third-order and in lllql where an elegant 
— »: , and non-perturbative non-linear anisotropy were got. The non-linear anisotropies have also been analyzed in lUU \12R base on 
f- >) ' the covariant approach to cosmological perturbations (see ll54ll for a recent review). Various secondary contributions to the non- 
Gaussianities have also been extensively studied, including the weak gravitational lensing and its correlation with ISW effect 
III7II21I - I26I1, wh ich is expected to be the dominant contamination of /^l^', correlation between lensing and Sunyaev-Zel'dovich 



(SZ) effect Il22[l . inhomogeneous recombination li!27--29il . small-scale dark matter clustering 1130113411 . A systematic treatment of 
^^ ■ the transfer function on all scales, which involves solving the full Boltzmann equations through the recombination phase and 
5-H ' then from the surface of last-scattering until today, has been performed in Oll - BSll at second-order. 



In this note, we calculate the CMB temperature anisotropy up to the third-order in primordial curvature perturbation Cwhich 
can be viewed as non-linear generalization of linear-order relation AT/T = — C/5- We follow the same strategy in ll8l- tlOl[l3l[l5ll . 
First, we calculate the gravitational redshift of a given photon from the emission surface to today's observer, which will give 
the observed anisotropy in terms of metric perturbations -^ = ^[$, ^, • • • ]. Then by using the conservation of curvature 
perturbation C in the large-scale limit, we determine the initial conditions for the metric perturbations in matter-dominated era, 
more precisely, the values $e = ^e[C] ^nd 'I'e = ^e[C] on the emission surface. Combining these two procedures will give the 
final non-linear mapping from ( to ^^ : -^ — -^ [(] . Together with previous results of primordial non-Gaussianities of ^ got in 
the literature, our formalism is ready to make prediction of the final observed CMB non-Gaussianity. Obviously, since we match 
the metric perturbations in matter era directly to those on the emission surface after last-scattering, our formalism includes only 
the gravitational redshift of photons, without considering the dynamics of photon-baryon plasma. Thus our result is valid only 
for large-scale anisotropies, which enter the horizon after decoupling. 
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' See e.g. I5g| for a recent short review of methods and techniques in calculating and analyzing primordial non-Gaussianities from inflationary models and 
references therein. 



This note is organized as follows. In the following section, we describe the temperature anisotropy induced by the gravitational 
perturbations from the emission surface to the observers. We give the general expression for the temperature anisotropy up to 
the third-order in metric perturbations. In the next section, by using the conserved curvature perturbation and taking the large- 
scale limit, we determine the initial conditions for metric perturbations in matter-dominated era. Then we give the non-linear 
mapping from primordial curvature perturbation to the temperature anisotropy, in the Sachs-Wolfe limit. Finally we give a short 
conclusion. 

II. FORMALISM 

After decoupling, the CMB photon density remains as Planck distribution, which is determined by a single parameter — the 
photon temperature. The temperature shift of a Planck distribution of photons is exactly proportional to the energy shift of any 
given photon, i.e. T/lo = const, if there is no collision, which implies 

^ = ^, (2.1) 

where Tf and Ti are the final and initial temperature of the Planck distribution respectively, w/ and tUi are the final and initial 
energy of a given photon respectively. ( 12.1b is exact, which implies that in order to get the change in the temperature, we need 
to evaluate the change in the energy of a (any) given photon, which in our question is nothing but its gravitational redshift. 
Thus, the question becomes to follow the geodesic equation of a given photon from the last-scattering surface to us, taking the 
inhomogeneous spacetime background into account. 
We work in the generalized Poisson gauge 

ds^ = a^g^^dx^'dx" = c? {-e^'^dif- + 2a,dridx' + e-^^"^ e^'' dx'dx^) , (2.2) 

where^ Ci i = 0, jij is symmetric and satisfies jij^i ^ ju ~ (thus det eJ^^ = 1), a is the scale factor The energy of a 
given photon with physical momentum P^ measured by an observer with 4-velocity u^^ = v^ /a (normahzed as a^gf^uU^u'^ = 

CO =. -a'g^.u^P'^ =. -g.^^u^p" , (2.3) 

where p^ is the momentum associated with the conformal metric g^^ (note p^ = P^/a^). Under the perturbed metric ( I2.2l i and 
using normalization for u'^ and p^^ to express w° and p" in terms of u' and p*, uj in terms of all relevant components takes the 
form; 



- ( \l 9i3P'P' \l 9iji'^'"3 + 1 - gjjW^J 



(2.4) 



with 



-2*„7,, , „-2* 



e7u +e-^*cr,cr,. (2.5) 



From (I2.4l l. on the background level u) = \p^\/a = p^ /a, which we will normalize by setting p'' = 1 in the following. The 
explicit expression for u^ (or v"^) requires details of the dynamics during recombination. 

Assuming the "intrinsic" photon temperature anisotropy at emission point x^ in the direction rXg takes the form 

T (^r]e,x],; nl) = T (^rje, x],; n*) e^v'='^<!'"<!J, where ije is the constant conformal time of emission, e.g. when the last-scattering 
takes place. The temperature measured by an observer at Xq and in the direction n is given by 

AT , ^ T{xl,n^)-T{xl,n^) a^A^^!?^^!?^ .(,.,.-„•) _ . ,. .. 

T T(xl,n^) aeU}e{ih,xl;p'e) 

where we have used To = ^T^ — —Te- From ( I2.6l l it is clear that the observed anisotropy comes from two aspects: the intrinsic 
anisotropy r on the emission surface which depends on the dynamics during recombination and the gravity theory, the other is 
the gravitational redshift-induced anisotropy from the emission surface to the observer ao'^o/{0'e^e), which is purely kinetic 
and is independent of the theory of gravitation. The intrinsic anisotropy r is highly model-dependent and a full treatment needs 
solving the set of Boltzmann equations up to the third-order which is beyond the scope of this note. On super-Horizon (the sound 
horizon on the last-scattering surface) scales where microscopic physics is irrelevant, a simple and non-perturbative expression 
has been got [16]: r ~ ~f ^^ which is adequate for our purpose. 

In the following, first we derive the gravitational redshift-induced anisotropy up to the third-order in terms of metric perturba- 
tions in ( I2.2I 1. then use the Einstein equation to determine the initial condition at matter-dominated era, i.e. the values of metric 
perturbations in terms of the conserved curvature perturbation Q on large scales. 



~ Here and in what follows, two repeated lower or upper spatial indices are contracted by Sij 

^ As is well-known, metric 3^,1/ = e^°'g^v and g^^ have the same null geodesies, but parameterized by different affine parameters A and A with relation 

d\ = e^"d\. Thus the energy in metric g^^ can be expressed as oi = —g/jiyuf'p'^ = —g^^^u^^p" , where p** = dx^ jdX and p^ = dx>^ jdX are momentum 

in metric g^^ and gp^y respectively. 



A. Perturbed photon energy 



After decoupling, the photons propagate freely and thus the physical content of its Boltzmann equation is completely encoded 
in the photon geodesic equation, which is much simpler to deal with. Null geodesies in perturbed spacetime has been investigated 
long before IIt UioI [o . [Isll . Although the whole thing we need to do is simply to follow the redshift of a photon, in a practical 
calculation, the complexities arise from several aspects. First, the energy to should be evaluated at the "real" emission point 
a;'(Ae) rather than at the "virtual image" at a;*(Ae) = (Ae — Xo)n'^- Here A is the affine parameter along the photon geodesies, 
Ae (Ao) is the corresponding values at emission surface (observer). Second, to determine the real position of emission we need 
to follow the photon geodesies, which we are able to solve only perturbatively. 

The expansion of frequency w around the background emission point x^ and the direction ?7* = — p* is straightforward ^- 
U^ll^uM- Here we simply report the corresponding expansions of ( I2.4l i up to the third-order in the metric perturbations ( I2.2l i. 
which involve the spatial components w* and p\ At the linear-order: 



"* ^ * + o7y"-*"^ 



where throughout this note we take the expansion of variable Q as Q ~ Q + Q(i) + Q(2) + Q(?,) + 

«^(2) = («(2) -P(2) + ^(l)P(l)) "' + 2*' + *"' (P(l) " Hi)) + 2 (^« " "W)' 



(2.7) 



At second-order: 
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^H,)^ 



(2.8) 



and 



auj 



(3) 



«(3) " P(3) + (^-(1) + 2^(1)"') [dj^h + 2 (^(1) + ^(1)"'') ^J^fe"(i) ) + (^' 

+ Ai)P\2) + (a^('2) - a^wJ'a)) {n^djv\^) +P{i)) - 2 (^w) ^(1) 



(2) 



(l)P(l) j djvl^) 



* - 'i^P(i) + i^lkin^n' I n* - 7,jnJ + p\^-^ - v\^-^ 



P{2} - ^hph 



+ («(i) - P(i) + l^ln' - 2*n') (vl^^ + [x^^ + x^i^n^") 9,«^i)) 



„0 „fe _0 „0 „k 



x-(2) + a;(2)"-^ - a:^i)P(i)"^ " a^(i)Pa)) ^fc + 3 (^W + ^(1)"'') (^(i) + ^(1)"') ^^^^' (* " 2^^^"^ 



i„3 



+ l^'PiD (P(i))' - I (">(!))' - '^ (^^"0' - ^*' - *' (^">(i) - ^7.,'^'"^ - 2i;^i)n» 



2i"P(i; 



2 l"(i: 



+ ^(4«(i)-P(i)JP(i) 



Y^ bv^'n^ 



1 r 
2 



.V^V 



2 (n'^>fi)) - 4 (pfi)) + 87fc;n''>|i) + (^-fkin'^n^ - Sn'^pfi)) (7mn?^™r^") - 27 
P(i) - 2n^pfi)n') p^;) + f|i) (i;^i) " 2^1;^)] 7,, + - (nVfi)^^-'' + 2vQ jf^n' + ^l%n'n^ 
+ (f(i) + 2:^(1)"^) { * - 2^*^"'"' + "' (P(i) - 2wii)) a,.* - -a^ (7yn*nJ)' 

(2.9) 
In (I2.7b -( |2!9l l, a dot denotes derivative with respect to A and subscripts "(j)" denote the orders in metric perturbations and all 
quantities are evaluated on the background emission point (virtual image)"*. In deriving the above results, we have used the 



"* )2.7t - fz!9t are essentially equal to (e.g.) eq.(2. 17)-(2. 19) in llal . The differences come from 1) here we have replaced v'^ and fe" in terms oft)', k^ and metric 
perturbations through the constraints v^ = —1 and k^ = 0, 2) we use <J> and 'J' rather than (p and -!/> which are defined as 1 + 2</> = e^*, 1 — 2ip = e~^* 
and 3) at this point we have not expanded $ = $(1) + ^12) + ^(3) s''-- 



expansion of Ag around its background value: Ag = Ag + A(x) + A(2) where A(i) = —x?i) (Ag) and A(2) = ^'■'^%) i^e) + 
^^1) (^'=) P'ii) i^^)- ^^^^ ^^^ be got by perturbing a::''(Ae) = S°(Ae) = rje, which is the definition of the emission surface as 
intersection of past light-cone of the observer and the spatial hypersurface at constant 1]^. 

B. Photon geodesies 

( I2.7b -( l2.9t can be fully determined when the perturbed photon geodesies is solved. Geodesic equation in the conformal metric 



p^ 



Kyp'^ 



0, 



(2.10) 



where p^ = dx'^ /d\, a dot denotes d/d\, Vt^ is the connection associated with the conformal metric g^^. 

Since i n (I2.7li - (l2.9l l we have expressed the photon energy in terms of spatial momentum p\ it is adequate to solve p* perturba- 
tively JTI-mX [l3l [isll . In order to make ( 12.10b a close set of equations for x*, we also need the expression for a;" in terms of x' . 
This can be done by perturbing the constraint p^p^ = around the background geodesies, which yields (up to the second-order 
in metric perturbation) 



P{1) 



P(l) 



A, 



(2.11) 



with 



A = $ + * 



Ctu' - -^ijn'n^ 



(2.12) 



and 



P(2) 



Pm 



1 



'>(i)P(i) 



<i)^M^+o(* + *) +2*'^'"' 



$ + '!' + --ijkn^v!" 



-($ + *) 7«j - ^(T^cro 



liji 



PW 



r'r,3 - - 



Jtjn'n^) 



(2.13) 



Integration of (12.11b and (12.13b with respect to A along the background geodesies a::^(A) will give x^^ and x%) ™ terms of xj-^-, 
and x%) respectively. 

Having deriving the general expressions for perturbed photon energy ( I2.7l )-( |2!9] l, in the following, we restrict ourselves to 
the large-scale limit. This is mainly because that in our formalism, to eventually determine the observed anisotropy in terms 
of conserved primordial curvature perturbation C,, we use the values of metric perturbations at matter-dominated era as initial 
conditions for gravitational redshift of photons rather than for the full Boltzmann equations. Thus, our results are only valid 
for the large-scale perturbation modes, which enter the horizon after decoupling and never affected by the microphysics. In the 
following we neglect the vector and tensor metric perturbations, not only because on large scales vector and tensor modes are 
subdominant but also calculation involving these modes up to the third-order is rather cumbersome. We also assume the the 
observer is comoving with the emission point, i.e w* = 0. Under these assumptions, a non-perturbative approach to the null 
geodesies has also been developed in ||15ll. 

Following the logic in I7 h10L[13L[15Ii . for spatially-flat FRW background, the set of perturbed geodesic equations are p^,. = 
a^(«) = /(„)' with 



/('l) = 


- -r,,(,)P''p% 




(2.14) 


/(*2) = 


= -r;.(2)PV - 2r;,(i)p^ifi) - xl^^dxT^.^^^p^p'' , 




(2.15) 


/(3) = 


~- ^^'pa(3)P''p'^ ^ '^^la{l)P''i'{2) ~ '^^la{2)P'' ^"{l) ~ ^'pa(l)^(l)^'{l) 








'^^iDdxri^iDP'ih - Uii)dxr;,i2)+^i2)dxrp,(i) + l^iiy^hdxdrri^^r) 


)?'?"■ 


(2.16) 



The perturbed Christofel symbol can be read from (IB7I ). After some manipulations, we can solve, at first-order in <!> and ^, 



P(i) 
Ph 



-2$ + /i, 



(2.17) 
(2.18) 



with 



dXd.A, 



K 



h 



dXA', 



where /! = $ + *. 



(2.19) 



Here and in the following we will frequently use the trick: (e.g.) ^ = (^' — n'9i^) 



At second-order iSHl 



with 



At the third-order [fisll . 



with 



^(2) = -2n' (a;fi)9^* + ^') - 2^11 + I^ 
2$2 - 2xf^.a,,$ - 2/1$ + II + /2, 



P(2) 



II = 



dA 



dA 



2($-/i)-xf,)a^ 



(a;^,)9^A' - 2$ A') 



ftA, 



^(3) 






d\ 



2 ($ - * - /i) xl'^^a^a,^ - x|'2)9^a,A - -xl^-^xl^^dxdrd.A 



2 ($2 _ 2*$) + "ill + 2/2 + 4 (* - $) h - 2a;{\.a^$ ) 9^^ 



(2.20) 
(2.21) 

(2.22) 
(2.23) 



(2.24) 



(2.25) 



by 



Finally, after plugging (I2.17b -( I2.24| | into ( I2.7b -( l279l l. the perturbed photon frequency up to the third-order in $ and ^ is given 



aaj(i) = -$ + /i, 



(2.26) 



and 



$ 



auj, 



(3) 



a,$ 



■^(2) 



$2/, 



1$-^ 



$/i + /2 - x^ijA' - (:r|i) + ^I'l^n^) 9,$, 



^ + ^^ + ll + h + h (2/2 - x?^') + $ x-^,) A' - Ii^ - /2 - .r^2)^' - ^ ^?i) ^" 



((a + * 



"^wj^w 






-x\2)+4i)Il--^\i){'^-h) 



In ( I2.26b -( l2728] l, all quantities are evaluated at the background emission point at xf'{\e) = (Ae — Aq)?!*. 



(2.27) 



(2.28) 



m. TEMPERATURE ANISOTROPY UP TO THE THIRD-ORDER 



A. Relation with primordial perturbations 

Having derived the perturbed photon energy in terms of the metric perturbations at emission surface, the next goal is to relate 
the temperature anisotropy to the primordial curvature perturbation C,, which encodes the information in the very early universe 
and is the most frequently used variable in evaluating the primordial non-Gaussianities in the literature. It is well-known that 
on large scales and for adiabatic perturbation, there is a non-pert urbative and conserved quantity which can be identified as 
non-linear curvature perturbation in uniform-density slices 145144811 . defined as^ 



T 1 r dp 

■ij-p 3(p + p) 



(3.1) 



Conserved and gauge-invariant C beyond the linear theory has also been constructed perturbatively in 11491 - 15 ID . The conservation 
of C will allow us to relate the primordial era when modes exit the horizon during inflation and the era when modes re-enter the 



^ The existence of non-perturbative conserved perturbation has also been derived by using covariant formalism I53 , l53ll . where the corresponding quantity is 
defined as a co-vector: C,a = daCc — "daP, a = o / drX'au'^ is the local expansion. See 15411 for a recent review. 



horizon, which is just the time of emission rjf, for our purpose. Our next task is to set the initial conditions for $ and ^P in the 
matter dominated era up to the third-order in the primordial curvature perturbation (. 
In matter-dominated era (p = 0), (13. lb can be integrated to give 



c = -* 



li Pn 

olll — 
3 Pn 



(3.2) 



On large scales, the matter density /?,„ can be related to metric perturbations through Einstein equation as (see Appendix (IB 11 1) - 

dni)): 



Pn 



-2* 



(3.3) 



which implies lITsllTq] 



c 



-* 



-$. 



(3.4) 



(I3.4l l is a non-perturbatively relation among C, and $, ^ on large scales during matter era. 

At linear order, for fluid without anisotropic stress, $ n;i = ^(i). which gives the well-known relation C = ^3*^(1)- However, 
<1> 7^ ^P at non-linear orders even for perfect fluid Il39t 14011 (see also ll4ll - l44ll for the discussion of evolution of higher-order 
cosmological perturbations). From the traceless part of {i — j)-component of Einstein equation and using (0 — 0) and (0 — i) 
components to express p and Ui in terms of metric perturbations (see Appendix (IBllb - (IB12b ). we are able to write a non- 
perturbative constraint between <!> and 'f on large scales during matter-dominated era [ 15, 16], 



d^{-^ -^)^^{d'^^f 






a$aa2$ _ 2a*a<92^ 



(3.5) 



In deriving ( 13.5b we have neglected higher-order spatial derivative terms since we are focusing on large scales. From ( 13.5b ^ can 
be solved up to third-order in $ as 



* = $ + 9" 



5 (92$)' + ^ (a,a,$)' + ^a,$9,92$ 



(3.6) 



This is the generalization of the linear-order relation ^ = $ up to the third-order in $. It is interesting to note the third-order 
part of (13.6b exactly vanishes. 

Combining ( 13.4b and (13.6b . it is now straightforward to solve $ = $[C] = $(i) + $(2) + ^(3) + 



$(1) = -3C, 



3{d\Y^ + id,d,cf + id,cd,d\ 



perturbatively to give 

(3.7) 
(3.8) 



and 



*(3) = - 



5^. 
125* 






(3.9) 



( I3.7b -( l379b give the large-scale initial condition for $ during the matter era, in terms of the conserved primordial curvature 
perturbation. In the above d~^ etc. can be understood in momentum space. From (I3.7I) -( [T9] | and (13.6b the corresponding initial 
condition ^P = '5 [C] up to the third-order in ( can also be easily get. 



B. Non-linear temperature anisotropy 



In the last part of this note, we will relate the observed temperature anisotropy ^ to the primordial curvature perturbation 
C, on large scales. To this end, we also need the intrinsic temperature at the emission surface in terms of metric perturbations. 
A fully treatment involves dynamics during recombination |31:r38ll . Here in this note, we take the large-scale non-perturbative 
expression found in II16I1 . where in matter dominated era with adiabatic assumption: Tg = Tge^**. Thus the large-scale 
temperature anisotropy up to the third-order in $ is given by ^ = {^-) ,-.-. + {^r') (2) '^ (t^) ni ^ ' ' ' ^^'■^ 



^ /(I) 
AT\ 



3 

l8 



ia,* (n'4) 



'(1) 



~3~ 



(3.10) 



4)^'' 



and 



/AT\ $3 $2j^ 



V T /f3> 162 18 



4)^' + ^*K)^'-^2)"^3 + ^4) 



a-^' ( 7i*x° 



^i)+2xl^A ~ 2hA' + x'l^^A'- 



a$ 



i (x^(2) + x°i)/i) + ^71* (34) + .tO^) ($ + 6A - 18a;°i) - 6/i)) + x^^ ($ - 3/i) 



(3.11) 



ia9j$ (x^d) + n'x^i)) (x^,) + n^x^i)) + 2n'd,<^' (x^^))^ 



In the following, in order to further relate ^21 to (, we take the SW contribution llTql where we neglect ISW and lensing 

6 3" — 1. Using (I3.7b -( ll!9] l. in momentum space, the non-linear mapping from ( to 

-|C' ^t the second-order 



contributions: ^ ^ f + fs + "^2 



■^ up to the third-order in C can be easily get. At linear order we find the familiar relation (■^) /, -, 



and third-order we find 
where the kernel 



73-^ if^-Pi-P2)l3[k;pi,P2)CpiCp2, 



{2^Y 
1 f d?pid?p2(Ppz 



3! 



{2nr 



S^{k-pi p2 -P3)l{k;pi,P2,P3)CpiCp2Cp 



Hk;PuP2) = -^ + ^ ^^fc^. 

7(fc;pi,P2,P3) = -j^ + [{l-giPi,P23))giP2,P3) + '2cyc\ic\ 



with pij = Pi + pj and 



9 {p, q) 



250 



l + 2P' + <-3-(^'"^')^ 



(p + qf (p + q)^ 



It is interesting to note that /3 — >■ 1/25 and g — ?> when p — >^ or q ^ 0, which implies that the non-linear m app i ng (13. 12I )-(|3 
would not contributes significantly in the so-called "squeezed" momenta configurations of non-Gaussianity lollsLllal . 
Using dASI l and ( IA6b . the bispectrum and trispectrum for ^ can be read as 



(3.12) 
(3.13) 

(3.14) 
(3.15) 

(3.16) 

13 



and 



B (fci, fc2, fcg) = a^Bf (fci, fc2, fcg) + {a^ Pq (fci) Pc (fc2) /? (fcs; fci, fca) + 2 cyclic) (3.17) 

T (fci^fc2, fca, ^4) =a^r^ (fei,fc2, ^3, ^4) 

+ a^B(^ {kik2, ki2) Pq (/ca) (3 (k^; ki2, fca) + llperms 

(3.18) 
+ a^Pc (^1) Pc. (fe) ^C (^3) 7 (fc4; -fei, -fe2, -fcs) + 3perms 

+ a^Pc (fci) ^C (^2) Pc (fcia) /3 (fcg; fci, fcig) /? (fc4; fc2, ^13) + llperms. 

Here P ^- an d T^ are primordial bispectrum and trispectrum for the curvature perturbation C respectively. 

In Ja Dal . the non-Gaussianities are conventionally characterized by non-linear relation of the Bardeen potential^ 

$ = $L + /nL * $! + gNL * $L, (3.19) 

where $l is the Gaussian part of $ and /nl and gNL are the so-called non-linear parameters, "*" denotes possible integration in 
momentum space^. To make contact with previous analysis and conventions in the literature, in the following we use the linear- 
order relations during matter-dominated era: -^ = — 1$ and if&L = |Cl, and make the ansatz for primordial non-Gaussianity: 
Q = Ql + |/nl * Cl + ^ffNL * Cl- After some manipulations, the non-linear parameters defined in ( |3.19t can be calculated as 



25 
/nl {k;pi,p2) = /nl ik;Pi,P2) - -ttP (k;pi,P2) , 



(3.20) 



* Here the Bai'deen potential $ should not be confused with the metric perturbation $ in )2.2t . Actually at linear order, <I>l = — $n j . 
' For example, /nl * *? = / "'^i^yr S^ (fc - Pi - P2) /nl (kipupi) ^L (pi) *L (P2). 



and 

50 r t 1 125 (3.21) 

- y [/3(fc;Pi,|fe-Pi|)/NL(l'^-Pih-P2,P3) + 2cyclicJ - — 7 (fe;Pi,P2,P3) ■ 

where /^^^ and g^^ ^^ non-linear parameters for primordial curvature perturbation C,, the functions /3 and 7 are given in ( 13.141 1 
and ( 13.15b . In ( I3.20l i- (I3.21| |. different contributions to the finally observed non-Gaussianity from primordial epoch and from 

AT 



non-linearity between C, and -^ are transparent. 



IV. CONCLUSION 



In this note, followin g th e approach developed in JSl- UOll . the general expression for the observed CMB anisotropy up to the 
third-order is calculated lUSll . In the Sachs-Wolfe limit, we derive the non-linear relation between the observed anisotropy and the 
conserved primordial curvature perturbation Q, up to the third-order in ^. ( I3.12b - (l3.13b can be viewed as non-linear generalization 
of familiar linear relation AT/T = — C/5. Our formalism is valid for large-scale anisotropies, which re-enter the horizon after 
decoupling. The results (I3.17b -( |3T8] | clearly show the different contributions to the observed non-Gaussianity from primordial 
non-Gaussianities in C, and non-linear mapping from C, to -^ due to gravitational perturbations. We also derive the non-linear 
parameters /nl and ^nl (eq.( l3.20l )-( l3.21l )). which enter the theoretical predictions for the angular bispectrum and trispectrum of 
CMB respectively. 

We do not expect the non-linear mapping in the SW limit ( I3.12b - (l3.13b would give a major contribution to the observed non- 
Gaussianities, especially comparing scenarios where large "primordial" non-Gaussianities can be generated (see IdSII for a recent 
review and references therein). However, the other secondary anisotropies, especially the correlation between lensing and ISW 
effect which we do not discuss in this note, is expected to give contribution to the final non-Gaussianity lll7ll2ll - l26ll . Another 
interesting issue is that, if an enhancement of non-linearity between -^ and C, is possible, like the enhancement of primordial 
non-Gaussianity by small Cs etc. in some inflationary scenarios. We wish to come back to these subjects in future investigations. 
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Appendix A: Non-Gaussianities from non-linear mapping 

Non-Gaussian variables can be get from non-linear mapping from Gaussian/non-Gaussian variables. The (JA^-formalism is a 
non-linear mapping from inflaton fluctuation 54) to the curvature perturbation C on super-Hubble scales. In general, the non-linear 
mapping in real space from a single variable Q io C, comes from three types: local products (e.g. JiV-formalism), products of 
local derivative terms, products of non-local terms. In all cases, the mapping can be written in fourier space as: 

1 r 1 r 

C,k = aQk + —^ j dpidp2Pk;pi.p2QpiQp2 + ^ / dpidp2dp3^k-p^^p^^p^Qp^Qp^Qp^ H , (Al) 

with dp = . ^3 , a is a fe-independent number and 

Pk-puP2 = (27r)^5(fc-pi -p2)/3(fc;Pi,P2), (A2) 

Ik;pi.p2.p3 = (27r)^5(fc-pi -P2 -P3)7(fe;Pi,P2,P3), (A3) 

where (3 {k;pi,p2) and J {k;pi,p2,P3) are normal functions. It is useful to note that f3 and 7 are symmetric with respect to all 
Pi's. 

It is more convenient to calculate the correlation functions of C = C ^ (C) since {() = 0. Straightforward calculation gives 
(subscript "c" denotes connected contribution): 



(A4) 



Ck^Ck2Cks) ^ {27:f S"" lY.kA B^{h,k2,h) , 

Ck,Ck2Ck,Ck,) ^ = {2TTf slY^kAr^ {kik2,k3, fc4) , 



with the leading contributions: 

Be (fci, k2, fca) = a^BQ (fci, fca, fcs) + {a'^Pg {ki) Pq (^2) P (fcg; fci, fca) + cyclic) , (A5) 



and 



T^ (fcl fc2, fcg, fc4) =a% (fel fc2, feg, fe4) 



a^JSg {kik2, ki2) Pq (fca) /3 (^4; ^12,^3) + Uperms 

a^Pg (fci) Pq (^2) Pq (fcs) 7 (^4; -fci, -^2, -fes) + Sperms 

a'^Pq (fci) Pq (fc2) Pq (^13) /3 (^3; fci, ^13) (i {kr, fc2, ^13) + Uperms 



(A6) 



with kij = \ki + kj\. In the above Pq is the power spectrum of Q, Bq and Tq are "intrinsic" bispectrum and trispectrum for Q 
respectively which vanish if Q is purely Gaussian. 

In 1I55I1 . based on the (5A^-formalism, Feynman-type diagrams were introduced to represent various contributions to non- 
Gaussianity from the non-linear mapping from 6(() to C, on large scales. In general this can be generalized straightforwardly for 
the non-Unear mapping jAll ). as we show in fig[T]and fig|2l 



I ^3 
t 



/3{fc3,fcl,fc2) 






Poih) Poih) 

FIG. 1: Diagrammatic representation of P{k\)P{k2)p{ks, k\, ^2) in ( jASt . 
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PQ(ki) 




A Pil'^i.ku.ks) 




X I(k4,-ki.-k2,-k3) 



ki 



I3{k3,ki,ki3) *-- 



PQiku) 



PcAkl) 



-« f-i{ki. k-i; kn) 



Poiks) 



FIG. 2: Diagrammatic representation of the second, third and fourth line in l lA6t . 



Appendix B: Perturbed metric and related quantities 

For the perturbed conformal metric 17^1, defined in ( 12.2b . the components of Christoffel symbol are 



-2* 

T°.=- 

^j 2Af2 






e '*'*' Gi 



-L n 



00 



2($+*) / 2* -7j„ > 



(Bl) 



(B2) 



(B3) 



(B4) 



i Oj - ^ jO 
2($+*) 



Ar2 



-e ^''=o-feO,<l> 



(e-^'" - ^ (e-^™a,„) (e'^-a,,)) [{-2^'e^^'' + (e^^-)') + e^* {d,ak - dk<J,)] 



(B5) 
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n-fc -^^^ [«'* (^.-^fc + ^^'^i) + 2*'e^^'= - (e^^'")'] 



1/,-.. - 



2* 



(B6) 
In the above N'^ = e^* + e'^'^e^^'-^aiaj, where N is the corresponding laspe function in ADM formalism. 

For metric ds^ = a^ (—e^'^drf + e~^^ dx^dx^\ the Christoffel connection is significantly simplified, with components: 



The corresponding components of Einstein tensor are 

Goo = 3 (H - *')' - e'^*+*^ ({d^-^f - 25^*) , 



(B7) 



(B8) 



Go. = 2 [a,*' + 9,$ (H - *')] , 



(B9) 



Gy = |e-2(*+*) [_2H' - (7^ - 2$' - 3^') [U - *') + 2*"] + {d,^f + d^ [^ - *)} J,, 

The above expressions are exact, which can be easily expanded to the desired orders. 
The (OO)-component of Einstein equation G^^, — T^^ gives. 



Pn 



3 (H - ^'f - e2(*+*) Ud.^f - 2^2*) g-2$ 



3^2 



(BIO) 



(BID 



where we used p,„ = ^i^- on the background level. From the (Oi)-component of Einstein equation we can solve 



1 + a ^e^ UiUi 



3H2-e2(*+*) {{d,^f ^28^-^ 



3-H2 _ e2(*+*) ( {d,^f - 202^ ) - 4H2e2(*+*) {d,^ 



(B12) 



where we have set ^' = 0. Finally, (IB lib and (IB 12b imply 

p2($+*) 



Pm — 2<I> 

Pm 



1- 



37^2 



{d,^Y -2d^^\ - 



4g2($+*) 



{d,<^Y 



9^2 



i-^i;^((a.*f-2a2* 



(B13) 



which gives, on large scales, Pml Pm — e 



-2* 
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